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ABSTRACT 

Infrared-dark high- mass clumps are among the most promising objects to study the 
initial conditions of the formation process of high-mass stars and rich stellar clus- 
ters. In this work, we have observed the (3-2) rotational transition of C 18 with 
the APEX telescope, and the (1,1) and (2,2) inversion transitions of NH3 with the 
Australia Telescope Compact Array in 21 infrared-dark clouds already mapped in the 
1.2 mm continuum, with the aim of measuring basic chemical and physical parameters 
such as the CO depletion factor (/d), the gas kinetic temperature and the gas mass. 
In particular, the C 18 (3-2) line allows us to derive /b in gas at densities higher 
(and hence potentially more depleted) than that traced by the (1-0) and (2-1) lines, 
typically used in previous works. We have detected NH3 and C 18 in all targets. The 
clumps have a median mass of ~ 244 M , are gravitationally bound, have an aver- 
age kinetic temperature of 17 K and possess mass, H2 column and surface densities 
consistent with being potentially the birthplace of high-mass stars. We have measured 
/d in between 5 and 78, with a mean value of 32 and a median of 29. These values 
are, to our knowledge, larger than the typical CO depletion factors measured towards 
infrared-dark clouds and high-mass dense cores, and are comparable to or larger than 
the values measured in low-mass pre-stellar cores close to the onset of the gravitational 
collapse. This result suggests that the earliest phases of the high-mass star and stellar 
cluster formation process are characterised by /b larger than in low-mass pre-stellar 
cores. On the other hand, /b does not seem to be correlated to any other physical 
parameter, except for a faint anti-correlation with the gas kinetic temperature. Thir- 
teen out of 21 clumps are undetected in the 24 /am Spitzer images, and have slightly 
lower kinetic temperatures, masses and H2 column densities with respect to the eight 
Spitzer-bright sources. This could indicate that the Spitzer-dark clumps are either less 
evolved or are going to form less massive objects. 

Key words: Molecular data - Stars: formation - radio lines: ISM - submillimetre: 
ISM - ISM: molecules 



1 INTRODUCTION 

An ever increasing number of observational evidences indi- 
cates that the earliest phases of massive star and stellar clus- 
ter formation occur within infrared dark clouds (IRDCs). 
These are dense molecular clouds seen as extinction fea- 
tures against the bright mid-infrare d Ga l actic background 
(e.g. Simon et al|200fl , Ragan et al. l2006l . l201ll . Rathborne 
et al. I2007L l20ld . Butler & Tan l2009h . IRDCs are charac- 
terised by very high gas column densities (10 23 — 10 25 cm -3 ) 
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and low temperatures 25 K), so that they are believed 
to be the place where most of the stars in our Galaxy are 
being formed. 

From an observational point of view, the spatial dis- 
tribution of the IRDCs in the Galaxy follows that of the 
molecular galactic component (with a concentration in the 
so-called 5 kpc molecular ring), and they are strong emitters 
of both far-IR/millimetre continuum and rotational molecu- 
lar transitions, especially those characterised by a high crit- 
ical density. In many of them, observations of the dense 
gas at sub-parsec linear scales revealed the presence of on- 
going star formation, both in isolated and clustered mode 
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(e.g. Beuther & Sridharan [20071 Zha ng et al. 20091, Fon tani 
et al. 120091 . Jimenez-Serra et al. l20ld . Pillai et al. l201lh . in- 
cluding clear signs of high-mass star formation like hot cores 
(e.g. Rathborne et al. 2008) and/or Ultracompact Hn re- 
gions (Battersby et al. l2010h . This demonstrates that IRDCs 
are indeed the birthplace of stars and stellar clusters of all 
masses. Therefore, the IRDCs in the earliest evolutionary 
stages are the best locations where to study the initial con- 
ditions of the star formation process and put constraints on 
current theories. 

Despite the identification of thousands of IRDCs, the 
number of studies devoted to unveiling their physical and 
chemical properties remains still limited. There are very 
few targeted studies, e specia lly in the southern hemisphere 
(e.g. Vasyunina et al. 2009, Miettinen et al. l201lh , where 
there are fewer groundbased facilities operating in the (sub- 
millimetre and centimetre domain than there are in the 
north. In particular, very little is known about the chemistry 
of IRDCs: studies suggest chemical compositions similar to 
those observed in low-mass pre-stellar cores (Vasyunina et 
al. [2009), including large abundances of deuterated sp ecie s 
(Pillai et al. 120071 : Pillai et al. l201ll : Fontani et al. l201lh . 
However, the amount of CO freeze-out, a key chemica l pa- 
rameter for pre-stellar cores (see e.g. Bergin & Tafalla 120071 
for a review) remains controversial: some w orks indicate high 
levels (factor of 5, Hernandez et al. 1201 lh of CO freeze-out, 
while others do not reveal significant CO depletion (e.g. Mi- 
ettinen et al. l201lh . 

In this work we present observations of the rotational 
transition (3-2) of the dense gas tracer C 18 0, performed 
with the Atacama Pathfinder Experiment (APEX) 12-m 
Telescope, and of the inversion transitions (1,1) and (2,2) 
of NH3 carried out with the Australia Telescope Compact 
Array (ATCA), towards 21 IRDCs with declination lower 
than —30° already mapped in the 1.2 mm continuum. The 
C 18 observations allow to compute the CO depletion fac- 
tor, while the ammonia inversion transitions can be used to 
derive t he tem perature in dense and cold gas (see e.g. Ho & 
Townes fl983h . Observations of N2H + and N2D + in some of 
the target sources, useful to derive the amount of deuterated 
fraction (by comparing the N2D + and N2H + column densi- 
ties), are also presented. In Sect. [2] we describe the criteria 
applied to select the targets. Sect.[3]gives an overview of the 
observations. The results are presented in Sect. 2] and dis- 
cussed in Sect. [5] Conclusions and a summary of the main 
findings are given in Sect. [6] 



2 TARGET SELECTION 

We selected 21 IRDCs from the 95 massive millimetre 
clumps detected by Beltran et al. (2006) in the 1.2 mm con- 
tinuum and non-MSX emitters (neither diffuse nor point- 
like). The targets were chosen according to these criteria: 
(i) source declination 5 ^ —30°; (ii) clumps isolated or hav- 
ing the emission peak separated by more than the SIMBA 
half power beam width to that of MSX-emitter objects, to 
limit confusion and select the most quiescent sources; (iii) 
clump masses above ~ 35M to deal with possible massive 
star formation. In this work we have recomputed the masses 
utilising the gas temperature derived from ammonia for each 
clump, and our results confirm that all targets are high-mass 



clumps (see Sect I4.6[) . The list of IRDCs is given in Table Q] 
The coordinates correspond to the peak of the 1.2 millime- 
tre continuum emission mapped by Beltran et al. (|20Q6h . 
In Table [1] we also give some basic information like the dis- 
tance to the Sun, the Galactocentric distance, and the (non- 
detection in the Spitzer-MIPS 24 /xm images. A comparison 
between the 1.2 mm continuum maps and the Spitzer 24 fim 
images of each target (except for 1 3039-6 108c6, for which 
the MIPS images are not available) is shown in Fig. IA1I of 
Appendix [Al 



3 OBSERVATIONS 

3.1 APEX 

Single-point spectra of the C 18 (3-2), N 2 H + (3-2) and 
N2D + (4-3) lines towards the sources listed in Table [1] were 
obtained with the APEX Telescope in service mode between 
the 20th and the 28th of June, 2008. The observations were 
performed in the wobbler- switching mode with a 150" az- 
imuthal throw and a chopping rate of 0.5 Hz. The receiver 
used for all lines was SHFI/APEX-2. The backend provided 
a total bandwidth of 1000 MHz. Details about the observed 
lines and the observational parameters are given in Table [2] 
The telescope pointing and focusing were checked regularly 
by continuum scans on planets and the corrections were ap- 
plied on-line. Calibration was done by the chopper-wheel 
technique. Spectra were obtained in antenna temperature 
units (corrected for atmospheric attenuation), TjJ, and then 
converted to main beam temperature units through the re- 
lation Tmb = jf^T% (F e ff and B e R are given in Col. 5 of 
Table [2]). The velocities adopted to centre the backends are 
given in Col. 6 of Table [1] For most sources we knew the ra- 
dial Local Standard of Rest velocities (Vlsr) from previous 
CS observations. The spectra of 16435-4515c3 (for which we 
did not have the CS data) were centred at Vlsr = 0; in any 
case the total bandwidth of the backend (~ 1000 km s _1 ) 
was larger than the velocity gradient across the Galaxy. 

3.2 ATCA 

We observed the inversion transitions (J,K) — (1,1) and 
(2,2) of ammonia at 23694.5 MHz and 23722.6 MHz (K-band 
at rsj 1.2 cm), respectively, with the ATCA0 towards all tar- 
gets in Table [1] The observations were performed between 
the 4th and the 8th of March, 2011, for a total telescope time 
of 48 hours. We used the configuration 750D, which pro- 
vides baselines between 31 and 4469 m. The primary beam 
was ~ 2.5' at the line frequencies. The flux density scale 
was established by observing the standard primary calibra- 
tor 1934-638, and the uncertainty is expected to be of the 
order of ~ 10%. Gain calibration was ensured by frequent 
observations of nearby compact quasars. The quasar 0537- 
441 was used for passband calibration. Pointing corrections 
were derived from nearby quasars and applied online. Atmo- 
spheric conditions were generally good (weather path noise 
~ 400 /im or better). 

1 The Australia Telescope Compact Array is part of the Australia 
Telescope which is funded by the Commonwealth of Australia for 
operation as a National Facility managed by CSIRO. 
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Table 1. Source list and detection summary. 



Source 


name a 


R.A. (J2000) 


Dec. (J2000) 


I 


b 


V LS R b 


d c 


D GC 


Mil 






h m s 


o 1 II 


o 


o 


km s _1 


kpc 


kpc 




08477- 


-4359cl 


08:49:35.13 


-44:11:59 


264.69 


-0.07 


8.5 


1.8 


8.9 


Y 


13039- 


-6108c6 


13:07:14.80 


-61:22:55 


305.18 


1.14 


-26.2 


2.4 


7.4 




13560- 


-6133c2 


13:59:33.04 


-61:49:13 


311.23 


-0.35 


-58.4 


5.6 


6.4 


N 


14183- 


-6050c3 


14:22:21.54 


-61:06:42 


314.03 


-0.52 


-42.6 


3.4 


6.6 


N 


15038- 


-5828cl 


15:07:32.52 


-58:40:33 


320.19 


-0.77 


-67.1 


5.0 


5.7 


N 


15278- 


-5620c2 


15:31:44.17 


-56:32:08 


324.07 


-0.73 


-49.4 


3.4 


6.1 


N 


15470- 


-5419cl 


15:51:28.24 


-54:31:42 


327.51 


-0.83 


-61.7 


4.1 


5.5 


ye 


15470- 


-5419c3 


15:51:01.62 


-54:26:46 


327.51 


-0.72 


-61.7 


4.1 


5.5 


Y 


15557- 


-5215c2 


15:59:36.20 


-52:22:58 


329.81 


0.03 


-67.6 


4.4 


5.2 


Y 


15557- 


-5215c3 


15:59:39.70 


-52:25:14 


329.80 


-0.00 


-67.6 


4.4 


5.2 


N* 


16061- 


-5048cl 


16:10:06.61 


-50:50:29 


332.06 


0.08 


-51.8 


3.6 


5.6 


Y 


16061- 


-5048c4 


16:10:06.61 


-50:57:09 


331.98 


0.00 


-51.8 


3.6 


5.6 




16093- 


-5128c2 


16:12:55.46 


-51:43:22 


331.77 


-0.86 


-97.3 


6.1 


4.3 


N 


16093- 


-5128c8 


16:12:49.63 


-51:36:34 


331.84 


-0.77 


-97.3 


6.1 


4.3 




16164- 


-4929c3 


16:20:24.51 


-49:35:34 


334.11 


-0.17 


-33.8 


2.6 


6.3 


N# 


16254- 


-4844cl 


16:29:00.89 


-48:50:31 


335.63 


-0.65 


-45.5 


3.4 


5.6 


Y 


16435- 


-4515c3 


16:47:33.13 


-45:22:51 


340.31 


-0.71 


0.0 


_h 


4.9 


N 


16482- 


-4443c2 


16:51:44.59 


-44:46:50 


341.24 


-0.90 


-43.3 


3.7 


5.1 


N/ 


16573- 


-4214c2 


17:00:33.38 


-42:25:18 


344.08 


-0.67 


-23.7 


2.6 


6.0 


Y 


17040- 


-3959cl 


17:07:58.78 


-40:02:24 


346.82 


-0.35 


-0.4 


16.5* 


8.5 


N?* 


17195- 


-3811c2 


17:23:00.30 


-38:14:58 


349.97 


-1.68 


-26.4 


3.6 


5.0 


Y 



a derived from the IRAS source name and the 'clump' number assigned by Beltran et al (120061 ): 
b from the CS (2-1) or (3-2) transitions (Fontani et al . l2005h : 

c kinematic distance from the Sun (Beltran et al 2006). In case of distance ambiguity, the near distance is adopted; 
d detection (Y) or non-detection (N) of the clump in the Spitzer-MIPS 24/im band; 
e point-like infrared source not at clump centre, probably associated with it; 

f point-like infrared source at the edge of the 3 a rms contour level of the millimetre continuum, probably not associated with it; 

9 diffuse infrared emission inside the 3 a rms contour level of the millimetre continuum, probably associated with a nearby infrared object; 

h the kinematic distance recomputed in Sect. 14.61 from the C ls O (3-2) line peak velocity is 3.1 kpc; 

* near distance smaller than 100 pc, so that the far distance is adopted; 

1 the millimetre map of the clump is incomplete because at the edge of the SIMBA field of view, thus one cannot conclusively assert if the 
clump is associated with an infrared source (see Fig.lATT). 



Table 2. Transitions observed with APEX and observational parameters. 



Transition 




HPBW 


AV b 


^eff/^eff 


T S ys 




GHz 


ii 


km s _1 




K 


C ls O (3-2) 


329.331 


18.9 


~ 0.22 


0.74/0.95 


200 - 300 c 


N 2 D+ (4-3) 


308.422 d 


20.2 


~ 0.23 


0.74/0.95 


160 - 300 


N 2 H+ (3-2) 


279.512 6 


22.3 


~ 0.26 


0.74/0.95 


130 - 160 



a line rest frequency; 
b spectral resolution; 

c for 13039-6108c6, 13560-6133c2 and 14183-6050c3, T sys - 600 K; 

d from the Cologne Database for Molecular Spectroscopy (CDMS), and it refers to the FiF = 55— ^44 
hyperfme component; 

e from Pagani et al. {2009), and it refers to the F±F = 4 5 ^ 3 4 hyperfme component, which has a relative 
intensity of 17.46%. 



The total time was broken up into series of 3-/5-minute 
snapshots in order to improve the coverage of the uv-plane 
for each target. As a consequence of this observing strategy, 
the integration time on source is variable, and generally in 
between ~ 30 mins and ~ 1 hour. The CABEq correlator 
provided two "zoom" bands of 64 MHz each, with a spec- 
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tral resolution in each zoom of 32 kHz (-0.4 km s" 1 at the 
frequencies of the lines). The ammonia lines were observed 
in one zoom band. The data were edited and calibrated fol- 
lowing standard tasks and procedures of the MIRIAD soft- 
ware package. After the editing and calibration in MIRIAD, 
the data were imported in AIPS. Imaging and deconvolu- 
tion were performed using the 'imagr' task, applying natu- 
ral weighting to the visibilities. The ammonia emission was 
detected only on the short baselines, thus we discarded all 
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baselines > 30 kA. In order to obtain images with the same 
angular resolution, and comparable to the APEX beam, we 
reconstructed the images with a circular beam of 20" of 
diameter for all sources, except for 17195-381 lc2 and 17040- 
3959cl, that have a poorer UV-coverage, and for which we 
get a beam of roughly 20" x 40". Moreover, 16254-4844cl 
and 16573-42 14c2 were observed only once, making the clean 
impossible and thus forcing us to use the dirty image decon- 
volved with the dirty beam to determine the spectrum. 

In this work the ATCA data will be used only to derive 
the gas temperature from the NH3 (1,1) and (2,2) spectra 
at the dust emission peak (in Sect. 14.5)) . A complete presen- 
tation and analysis of the ATCA observations will be given 
in a forthcoming paper (Giannetti et al., in prep). 



4 RESULTS 

4.1 Detection summary 

The detection summary is given in TableU C 18 (3-2) was 
observed and detected in all sources. Twelve targets were 
also observed and detected in N 2 H + (3-2). Ten clumps were 
observed in N 2 D + (4-3), and only clump 15557-5215c2 was 
marginally detected. This is, to our knowledge, the first de- 
tection of this line towards an IRDC. All targets were ob- 
served and detected in the ammonia inversion transitions, 
with the exception of 14183-6050c3 which was undetected 
in the (2,2) line. 



4.2 Line profiles and parameters 

All spectra of the C 18 and N2H + lines are shown in Ap- 
pendix E Fig. flTTI shows the N 2 H+ (3-2) and C 18 (3-2) 
transitions in the twelve clumps observed and detected in 
both lines. Fig. IB 2 1 shows the nine sources observed and de- 
tected in C 18 (3-2) only. 

Most of the C 18 lines are well fitted by single Gaus- 
sians. The results of these fits are reported in Table [4l 
Slightly asymmetric profiles deviating from the Gaussian 
shape can be noted in 15470-5419cl, 16061-5048cl, 16482- 
4443c2, 16573-4214c2 and 17040-3959cl. This could be 
due either to the superposition of multiple blended ve- 
locity components, or to high optical depth effects. Non- 
Gaussian wings in one (or both) side(s) of the line indicat- 
ing presence of outflows are noticeable towards two sources, 
08477-4359cl and 17195-3811c2 (see Fig. [Q). Four sources 
show multiple lines well-separated in velocity, 14183-6050c3, 
16093-5128c2, 16093-5128c8, 16435-4515c3. Be cause we 
had the previous CS observations (Fontani et al. 120051 ) as 
reference, in Table 3] we have identified the line associated 
with the star- forming region of our interest. The other com- 
ponents are almost certainly arising from clouds along the 
line of sight but not associated with the star-forming re- 
gion, given the large separation in velocity (from ~ 20 to 
-40 km s" 1 , see Figs. IBT1 and lB2l). 

The N2H + (3-2) lines roughly peak at the same position 
as the C 18 (3-2) ones and have comparable line widths. 
Because of the hyperfine structure present in the transition, 
it is not straightforward to derive conclusions on blended 
multiple velocity components. However, we note hints of 
possible secondary velocity components in all the lines that 
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Figure 1. C ls O (3-2) spectra of 08477-4359cl and 17195- 
381 lc2. The Gaussian fits are superimposed on the spectra 
and highlight significant emission in non-Gaussian wings. 



Table 4. C ls O line parameters derived from Gaussian fits. For sources 
with more than one velocity component, the line labelled as '-a' is the 
one associated with the IRDC (based on previous CS observations, 
Fontani et al. 120051) . 



source 






/ T M Bdv 


^LSR 


FWHM 


Tpeak 








K km s- 1 


km s _1 


km s _1 


K 


08477- 


-4359cl 




6.10(0.09) 


9.53(0.02) 


2.44(0.04) 


2.3 


13039- 


-6108c6 




1.4(0.1) 


-26.04(0.06) 


1.8(0.1) 


0.76 


13560- 


-6133c2 




1.6(0.15) 


-57.47(0.08) 


1.9(0.2) 


0.8 


14183- 


-6050c3- 


a 


0.9(0.1) 


-45.52(0.08) 


1.5(0.2) 


0.56 


14183- 


-6050c3- 


b 


0.57(0.06) 


-26.70(0.02) 


0.48(0.05) 


1.1 


15038- 


-5828cl 




0.9(0.1) 


-69.3(0.3) 


5.2(0.7) 


0.16 


15278- 


-5620c2 




2.5(0.1) 


-49.42(0.04) 


2.0(0.1) 


1.18 


15470-5419cl 




2.3(0.1) 


-58.74(0.07) 


2.45(0.2) 


0.88 


15470- 


-5419c3 




2.6(0.1) 


-61.38(0.06) 


2.9(0.1) 


0.84 


15557- 


-5215c2 




4.2(0.1) 


-67.57(0.03) 


2.28(0.08) 


1.7 


15557- 


-5215c3 




1.7(0.1) 


-68.47(0.08) 


2.4(0.2) 


0.67 


16061- 


-5048cl 




7.61(0.02) 


-67.28(0.03) 


3.58(0.06) 


2.0 


16061- 


-5048c4 




3.2(0.1) 


-51.99(0.03) 


1.94(0.07) 


1.6 


16093- 


-5128c2- 


a 


3.6(0.1) 


-96.19(0.05) 


3.4(0.1) 


1.0 


16093- 


-5128c2- 


b 


1.97(0.07) 


-63.46(0.03) 


1.57(0.07) 


1.2 


16093- 


-5128c8- 


a 


0.96(0.06) 


-96.85(0.08) 


2.3(0.2) 


0.39 


16093- 


-5128c8- 


b 


0.19(0.04) 


-68.9(0.1) 


1.0(0.2) 


0.18 


16164- 


-4929c3 




2.95(0.07) 


-33.80(0.02) 


1.76(0.05) 


1.6 


16254- 


-4844cl 




2.72(0.09) 


-44.98(0.03) 


2.20(0.09) 


1.2 


16435- 


-4515c3- 


a 


1.7(0.1) 


-35.57(0.04) 


1.6(0.1) 


1.0 


16435- 


-4515c3- 


b 


0.88(0.09) 


-52.09(0.06) 


1.3(0.2) 


0.62 


16435- 


-4515c3-c 


0.5(0.1) 


-24.4(0.1) 


1.5(0.4) 


0.31 


16482- 


-4443c2 




4.6(0.1) 


-43.35(0.05) 


3.3(0.1) 


1.3 


16573- 


-4214c2 




8.1(0.2) 


-27.29(0.03) 


2.76(0.06) 


2.8 


17040- 


-3959cl 




6.45(0.08) 


-15.87(0.02) 


2.40(0.03) 


2.5 


17195- 


-3811c2 




2.54(0.07) 


-26.41(0.02) 


1.47(0.06) 


1.6 



show similar features in the C 18 (3-2) line (15470-5419cl, 
16061-5048cl, 16573-4214c2) except for 16482-4443c2. 

In Table [5] we give the N2H + (3-2) line parameters: in 
Cols. 3 - 7 we list integrated intensity (J TMBdv), peak ve- 
locity (Vlsr), FWHM, opacity (r), and excitation tempera- 
ture (T ex ) of the N2H + (3-2) line, respectively. Because the 
rotational transitions of N2H + possess hyperfine structure, 
we fitted the lines using the METHOD HFS of the CLASS 
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Table 3. Detection 


summary: 


Y = detected, N 


= not detected. 


, - = not observed; 


Source 




APEX 




ATCA (NH 3 ) 
NH 3 (1,1) / (2,2) 


C 


18 O (3-2) 


N 2 H+ (3-2) 


N 2 D+ (4-3) 


no/177 /|orn„i 

U8477— 4o59cl 


Y 


Y 


N 


Y / Y 


13U39-ol08co 


Y 


Y 


N 


Y / Y 


135oU-ol33cz 


Y 






Y / Y 


14183-o050c3 


Y 






Y / N 


15L)c58-58z8cl 


Y 




N 


Y / Y 


15z78-5dzUcz 


Y 




N 


Y / Y 


1547U— 5419cl 


Y 


Y 


N 


Y / Y 


15470-5419c3 


Y 


Y 


N 


Y / Y 


15557-5zl5cz 


Y 


Y 


ya 


Y / Y 


15557-5zl5c3 


Y 


Y 


N 


Y / Y 


16061-5048cl 


Y 


Y 


N 


Y / Y 


16061-5048c4 


Y 


Y 


N 


Y / Y 


16093-5128c2 


Y 


Y 




Y / Y 


16093-5128c8 


Y 






Y / Y 


16164-4929c3 


Y 






Y / Y 


16254-4844cl 


Y 






Y / Y 


16435-4515c3 


Y 


Y 




Y / Y 


16482-4443c2 


Y 


Y 




Y / Y 


16573-4214c2 


Y 


Y 




Y / Y 


17040-3959cl 


Y 






Y / Y 


17195-3811c2 


Y 






Y / Y 



a marginal detection (see Sect. 14.4]) . 



Table 5. N 2 H+ line parameters and column densities. For optically thin transitions without well-constrained opacity 
(r = 0.1), to compute Af(N 2 H + ) we adopted a T ex = 6 K (see text). For sources with more than one velocity component, 
the line labelled as '-a' is the one having the peak velocity consistent with that of the clump of our interest, as in Tabled 



source 


vel. range 


J T M Bdv 


^LSR 


FWHM 


T 




N(N 2 H+) 
(xl0 13 )cm- 2 




km s _1 


K km s- 1 


km s _1 


km s _1 




K 


08477-4359cl 


5.4,14.8 


6.02(0.07) 


10.53(0.02) 


2.87(0.05) 


0.1 


6 


2.93(0.03) 


13039-6108c6 


-30.0,-23.6 


0.66(0.08) 


-26.2(0.1) 


1.6(0.4) 


0.1 


6 


0.32(0.04) 


15470-5419cl 


-66.1,-54.8 


5.2(0.1) 


-58.82(0.03) 


2.4(0.2) 


12(2) 


4.8(0.3) 


54(1) 


15470-5419c3 


-67.1,-56.8 


9.4(0.1) 


-61.38(0.03) 


3.14(0.02) 


0.1 


6 


4.6(0.05) 


15557-5215c2 


-74.1,-62.6 


17.2(0.2) 


-66.77(0.02) 


2.01(0.05) 


12.7(0.8) 


8.0(0.5) 


27(2) 


15557-5215c3 


-71.8,-66.2 


1.12(0.07) 


-68.69(0.06) 


1.61(0.09) 


1.4(0.2) 


4.6(0.2) 


1.8(0.1) 


16061-5048cl 


-72.6,-59.5 


25.3(0.2) 


-66.20(0.02) 


3.14(0.09) 


9.1(0.8) 


9(1) 


29(3) 


16061-5048c4 


-56.3,-48.3 


2.13(0.08) 


-51.98(0.05) 


1.29(0.09) 


7(1) 


4.5(0.3) 


18(4) 


16093-5128c2-a 


-102.1,-92.1 


1.28(0.06) 


-97.20(0.08) 


3.1(0.2) 


0.1 


6 


0.62(0.03) 


16093-5128c2-b 


-68.0,-60.2 


0.63(0.06) 


-63.21(0.08) 


1.7(0.3) 


0.3(1.5) 


6 


0.38(0.04) 


16435-4515c3-a 


-53.0,-49.8 


0.08(0.02) 


-51.1(0.3) 


1.6(0.4) 


0.1 


6 


0.04(0.01) 


16435-4515c3-b 


-37.1,-32.5 


0.23(0.02) 


-35.2(0.1) 


1.5(0.3) 


1(0.2) 


3.5(0.7) 


2.7(0.2) 


16482-4443c2 


-47.0,-38.6 


3.68(0.07) 


-42.57(0.03) 


3.02(0.09) 


0.1 


6 


1.79(0.03) 


16573-4214c2 


-32.5,-22.2 


9.52(0.08) 


-27.06(0.01) 


3.79(0.02) 


0.1 


6 


4.64(0.04) 



hyperfine components, most of which are blended due to the 
fact that the line widths are larger (~ 1 — 3 km s _1 ) than 
the separation in velocity of the components. However, the 
fit residuals are generally low indicating that the procedure 
provides good results despite the blending. 

4.3 Column densities of N 2 H + 

The N 2 H + column densities were derived follow ing th e 
method outlined in the Appendix A of Caselli et al. (|2002al ). 
Specifically, we adopted equations (A-l) and (A-4) for opti- 
cally thick and optically thin lines, respectively. The method 
assumes a constant excitation temperature, T ex - An estimate 



packagjfl- This method fits all the hyperfine components si- 
multaneously assuming that they have the same excitation 
temperature and width, that the opacity has a Gaussian de- 
pendence on frequency and fixing the separation of the com- 
ponents to the laboratory value. The line parameters listed 
in Table [5] have been derived from this method, except the 
integrated intensity that has been computed by simple in- 
tegration over the velocity range given in Col. 2. The line 
opacity is deduced from the intensity ratio of the different 

3 The CLASS program is part of the GILDAS software, devel- 
oped at the IRAM and the Observatoire de Grenoble, and is avail- 
able at |http://www.iram.fr/IRAMFR/GILDAS| 
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of T ex can be derived from the fitting procedure to the hyper- 
fine structure Q, but not for optically thin lines or lines with 
opacity not well- constrained. For these, we have assumed 
the average Tex (~ 6 K; Table [5} derived for the sources 
with well-constrained opacity. 

For all sources but 16573-4214c2 and 16482-4443c2, 
the clump diameter derived from the 1.2 mm continuum 
is comparable t o or l arger than the beam size (see Table 2 
in Beltran et al. 2006, see also Table [8]), so that we assumed 
unity filling factor. This factor was applied also to the un- 
resolved source 16482-4443c2. We stress that this is an ap- 
proximation, because the sources could be clumpy. However, 
due to the lack of observations at higher- angular resolution, 
neither in the lines observed in this work nor in other high- 
density gas tracers, the size of the effective emitting region 
cannot be determined. For 16573-4214c2, for which 6 S ~ 7" , 
we applied a correction factor l/r] v — (0% + HPBW 2 )/0q ~ 
10. 

The N2H + total column densities are given in Col. 8 of 
Table [5] and span a range from 10 12 to 10 14 cm . These 
values are in good agreement with those derive d from the 
same line i n other IRDCs (Ragan et al. 120061 Miettinen 
et al. 1201 lh , as well as in ot her m assive starl ess an d star- 
forming cores (Fontani et al. 1201 lL Chen et al. 120111 ). 

4.4 Upper limits on the column density of N2D + 
and on N2D + /N2H + column density ratio 

The N2D + (4-3) transition was marginally detected at ~ 
3.2cr rms only towards 15557-5215c2. This represents, to 
our knowledge, the first detection of this line in an IRDC. 
The spectrum is shown in Fig. [2] For this source we have 
computed the N2D + column density fitting the line with a 
Gaussian (the results are given in Table [6]) , and t hen fo llow- 
ing the approach of Appendix A in Caselli et al. (2002a|) to 
derive the N2D + column density for the optically thin case. 
This line is also a blend of hyperfine components, therefore 
the line width derived from the Gaussian fit is an upper 
limit to the intrinsic value. However, because of the poor 
signal-to-noise ratio of the spectrum, METHOD HFS did 
not provide reliable results. As T ex , we used that derived 
from N 2 H + (3-2) for this source (see Table 0. We find 
Af(N2D + )^ 8 x 10 11 cm -2 and the corresponding deuterated 
fraction (column density ratio Ar(N 2 D + )/Ar(N 2 H + ) = Arac) 
is rsj 0.003 (Table [6]). This value is consisten t with studies 
performed in both IRDCs (e.g. Miettinen et al. 2011, Chen et 
al. 120111 ) and massive star- forming clumps co ntaining more 
evolved objects (e.g. Fontani et al. 120061 1201 lh . 

For all the other undetected sources, we have derived 
upper limits on the deuterated fraction in between 0.003 and 
0.5. These upper limits were computed from the integrated 
intensity upper limits assuming the lines to be Gaussian 
from the formula f TyiBdv = — ^ v TF A e ^ k . We adopted the 

J 2A/ln2/7r Mt5 

3cr rms level in the spectrum as peak temperature T^^ k , and 
the line width measured in the detected source 1 5557-52 15c2 
(Table E} as AV. The £> frac upper limits are comparable to 
the values commonly measured in massive clumps (Fontani 

4 See the CLASS user manual for the deriva- 
tion of T ex from the output parameters: 
|http://iram.fr/IRAMFR/GILDAS/doc/html/class-html/class.html/| 
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Figure 2. N2D+ (4-3) spectrum of the only source detected 
in this transition, 15557-5215c2. 



et al. 2006, Chen et al. l201ll . Miettinen et al. 12011 ). so that 
our sensitivity does not allow to conclude if our targets have 
deuterated fraction similar to or lower than the other clumps 
in the literature. We believe that the very low detection rate 
in the N2D + (4-3) line is likely due to the high densities re- 
quired to excite this transition: having a critical density of 
~ 3 x 10 7 cm -3 , this line is expected to come from very com- 
pact regions, thus suffering enormous beam dilution effects. 

4.5 Rotation temperatures and NH3 column 
densities at dust emission peak position 

We have extracted spectra of the NH3 (1,1) and (2,2) tran- 
sitions from the ATCA channel maps towards the positions 
given in Table [TJ The spectra obtained this way were then 
imported in CLASS, and fitted using METHOD NH 3 to fit 
the (1,1) lines, which takes into account the line hyperfine 
structure, similarly to METHOD HFS (see Sect. 02]). This 
was also used for the (2,2) emission, even though the hyper- 
fine components were not always visible, to obtain an upper 
limit for r. 

In Table we list the rotation temperatures (T ro t), ki- 
netic temperatures (Tk), and total column densities of am- 
monia (A^nh 3 ) derived from the NH3 (2,2)/(l,l) line ra- 
tios, as well as the line peak velocities. We have derived 
T r ot and A^nh 3 from the output parameters given by the fit- 
ting method outlined above, and usi ng th e formulae given 
in the Appendix of Busquet et al. (|2009h . The formulae, 
which result from the discussion in Ho & Townes (1983, Eq. 
(4)), have been derived assuming that the transitions be- 
tween the metastable inversion doublets are approximated 
as a two-level system, and that the excitation temperature 
and line width are the same for both NH3 (1, 1) and NH3 
(2, 2). Note that the assumption of a two-level system is 
reasonable because transitions between the metastable in- 
version doublets are usually mu ch fas ter than those of other 
rotational states (Ho & Townes 1983). Tk was extrapolated 
from T r ot following the empirical approximation method de- 
scribe d in T afalla et al. (J200J, see also Eq. (A. 5) in Busquet 
et al. I2009T ). We find kinetic temperatures in between 13 
and 25 K, with both mean and median Tk of ~ 17 K. These 
values are consistent, within the errors, with t ypical tem- 
peratures measured towards IRDCs (Pillai et al. I2006T ) from 
ammonia. 



CO depletion in infrared dark clouds 

Table 6. N 2 D+ (4-3) line parameters and N 2 D+ column density in 15557-5215c2 
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/ T MB dv 


^LSR 


FWHM 


^"peak 


iV(N 2 D+) 


7V(N 2 D+)/7V(N 2 H+) 


K km s -1 


km s _1 


km s _1 


K 


(xlO n )cm- 2 




0.14( 0.03) 


-68.1(0.2) 


2.0(0.5) 


0.07 


8(2) 


0.003(0.001) 



Table 7. Rotation temperature, kinetic temperature, total 
NH3 column densities and peak velocities derived from the 
NH3 inversion transitions observed with the ATCA. The er- 
rors on Trot and Tk are given in parentheses and are computed 
from the propagation of errors. The uncertainty on the column 
densities includes the error on the flux calibration, and is es- 
timated to be of the order of 20 - 30%. 



Source 


Trot 


T k 


%H 3 


^peak 


K 


K 


10 15 cm" 2 


km s _1 


08477- 


-4359cl 


16(2) 


19(3) 


1.06 


+11.18(0.02) 


13560- 


-6133c2 


15(9) 


17(10) 


4.61 


-56.32(0.05) 


14183- 


-6050c3 b 








-44.36(0.04) 


15038- 


-5828cl 


14(3) 


15(4) 


2.97 


-68.72(0.04) 


15278- 


-5620c2 


17(2) 


20(3) 


3.52 


-48.69(0.02) 


15470- 


-5419cl-a 


16(2) 


18(3) 


4.23 


-60.72(0.03) 


15470- 


-5419cl-b 


13(9) 


14(9) 


3.20 


-57.88(0.06) 


15470- 


-5419c3 


17(1) 


19(2) 


3.61 


-61.02(0.01) 


15557- 


-5215c2 


19(2) 


23(4) 


4.62 


-67.11(0.01) 


15557- 


-5215c3 


14(3) 


15(3) 


4.47 


-68.53(0.01) 


16061- 


-5048cl 


20(3) 


25(5) 


4.44 


-66.76(0.01) 


16061- 


-5048c4 


12(2) 


13(2) 


5.28 


-51.62(0.01) 


16093- 


-5128c2 


14(5) 


16(5) 


3.87 


-97.09(0.05) 


16093- 


-5128c8 


14(6) 


15(7) 


1.09 


-95.80(0.04) 


16164- 


-4929c3 


13(6) 


14(6) 


2.59 


-32.80(0.01) 


16254- 


-4844cl c 


18(4) 


22(6) 


4.77 


-38.45(0.05) 


16435- 


-4515c3 


11(4) 


12(5) 


1.84 


-51.69(0.04) 


16482- 


-4443c2 


15(3) 


16(4) 


0.78 


-42.35(0.03) 


16573- 


-4214c2 c 


15(4) 


17(5) 


3.32 


-25.09(0.04) 


17040- 


-3959cl 


16(3) 


18(4) 


6.81 


-15.62(0.02) 


17195- 


-3811c2 


15(1) 


17(2) 


1.13 


-25.51(0.01) 



a derived from the (1,1) transition. The uncertainty (between 

parentheses) is that given by the fitting procedure; 

b undetected in the (2,2) transition (see Table O; 

c observed only for one cycle with bad UV coverage. The 

values listed have been derived from the spectra extracted 

from the dirty map deconvolved with the dirty beam. 



4.6 



H 2 masses and other physical parameters 
from 1.2 mm continuum emission 



Gas masses were calculated by Beltran et al. (|2006h for all 
clumps from the 1.2 mm continuum integrated flux den- 
sity, assuming optically thin emission and a reasonable dust 
temperature of 30 K for all clumps. This latter assumption 
was due to the fact that a temperature estimate for each 
clump was lacking. In this work we can profit from the tem- 
peratures derived from ammonia (Col. 3 of Table [7|) and 
recompute the H 2 masses assuming that the dust temper- 
ature equals the kinetic temperature. This method implies 
coupling between gas and dust, which is a realistic assump- 
tion for gas densities as those of our clumps (^ 10 4 cm -3 ). 

The gas mass, M, has been derived using Eq. (1) in 
Beltran et al. (|2006l ). and adopting the same assumptions, 



except the dust temperature for which we have utilised Tk 
derived from ammonia. From M, for each source we have 
calculated the source averaged H 2 volume and column den- 
sities, riH 2 and 7Vt(H 2 ), assuming the clumps to be spherical 
and homogeneous. The H 2 column densities, iV p (H 2 ), have 
been estimated also from the 1.2 mm contin uum peak flux 
using the equations in Beuther et al. (2002) and adopting 
the same assumptions made to derive the mass (same /?, dust 
mass opacity and gas-to-dust ratio). These represent aver- 
age values in the telescope beam of the 1.2 mm continuum 
observations (~ 24"), and we will use these estimates to de- 
rive the CO depletion in Sect.EO because the APEX C 18 
observations have a comparable angular resolution (~ 19"). 
Also, we have computed the source averaged surface den- 
sity, S(H 2 ) = 4M/7r(6 s d) 2 , where d and S are the source 
distance (see Table [T]) and angular diameter (see Table [8]), 
respectively. For 1 6435-45 15c3, for which we did not have a 
kinematic source distance, we computed it from the velocity 
of the C 18 (3-2) line of the '-a' component (see Table 0J 
following the method explained in Fontani et al. (|2005l ). The 
metho d, which assumes the rotation curve of Brand & Blitz 
()1993l ). is valid for distances from the Galactic centre be- 
tween 2 and 25 kpc, and provides two kinematic distances, 
3.1 and 12.6 kpc. As for the other sources with distance 
ambiguity, we adopted the 'near' value, i.e. 3.1 kpc. 

The results are listed in Table [H where we give M, 
riH 2 , N(R 2 ) (both iV t (H 2 ) and iV P (H 2 )) and E(H 2 ). For com- 
pleteness, in Cols. 2 and 3 of Table [8] we list the full width 
half maximum diameters, S (deconvolved with the telescope 
beam of 24") adopted to derive riH 2 , 7Vt(H 2 ) and H(H 2 ), 
as well as the integrated flux densities used to calculate 
the masses, S u - Both parameters are taken from Beltran et 
al. ((2006). The mean mass turns out to be 439 M , and the 
median is 244 M©. As expected, these values are systemati- 
cally higher than those derived by Beltran et al. (|2006l ) given 
that the kinetic temperatures from ammonia (see Table [7} 
are systematically lower tha n th e representative value (30 K) 
assumed by Beltran et al. (2006). The clump- averaged col- 
umn densities are of the order of 10 22 - 10 23 cm , with a 
mean value of 1.6 x 10 23 cm -2 and a median of 8.5 x 10 22 
cm - . The volumn densities are of the order of 10 4 - 10 5 
cm -3 , with a mean value of 3.1 x 10 5 cm -3 and a median 
of 3.9 x 10 4 cm -3 . Finally, the mean surface density is 0.37 
g cm -2 , and its median value is 0.19 g cm -2 . The difference 
between mean and median values for all parameters is due 
to the clump 16573-4214c2, which has an angular diameter 
(~ 7") more than twice smaller than any other source of the 
sample, despite the comparable mass. 

The clump- averaged column densities A^(H 2 ) are gen- 
erally smaller than the theoretical threshold given by 
Krumholz & McKee (2008) to avoid cloud fragmentation, 
which is rsj 1 g cm -2 , corresponding to ^ 3 x 10 23 cm -2 . 
Also, the distance independent parameter S(H 2 ) is smaller, 
on average, than the theoretical values predicted for clumps 
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that are going to form high-mass stars or super star clus- 
ters, which are expecte d to b e larger than ~ 0.7 g cm -2 
(Chakrabarti & McKee I2OO5I . Krumholz & McKee 2008). 
Rather, the values measured in this work are consistent 
to those predicted for clumps that are going to form 
intermediate-mass stars and stellar clusters (Chakrabarti & 
McKee [jpoj ). We stress, however, that E(H 2 ) and i\T t (H 2 ) 
represent average values across the whole clumps, which 
could in reality be fragmented in smaller and denser cores. 
Therefore, our E(H2) and Nt(H.2) are to be considered as 
lower limits for the individual embedded cores. On the other 
hand, the column densities calculated from the 1.2 mm 
continuum peak flux, iVp (H2), are closer to the theoretical 
threshold proposed by Krumholz & McKee (2008), indicat- 
ing that the central regions of the clumps could more easily 
form massive stars. Moreover, we have compared the clump 
masses to the threshold proposed by Kauffmann & Fil- 
ial (2010) based on an empirical mass-radius relation which 
predicts that a cloud must exceed the relation M > M t h r = 
87OM (i?/pc) 1 - 33 to form massive stars. The mass thresh- 
olds M t hr are listed in Col. 9 of Table [8] As R, we have used 
half of the diameters listed in the same Table. We can see 
that our clump masses are all above the predicted threshold, 
with four exceptions: 13039-6108c6, 15557-5215c3, 16164- 
4929c3 and 16482-4443c2, for which, in any case, the mea- 
sured mass and the corresponding threshold value are sim- 
ilar, and for 16482-4443c2 the mass threshold is even an 
upper limit. Based on these results, we claim that the tar- 
gets have the potential to form massive stars. 



4.7 C ls O column densities and CO depletion 
factors 

The C 18 column density has been derived from line in- 
tensities assuming optically thin lines and LTE conditions. 
Under these assumptions one can demonstrate that the total 
column density of C 18 0, N c i8 Q , is related to the integrated 
intensity of a ro tation al transition J — »• J — 1 according to 
(e.g. Pillai et al. 120071 ): 

N c iso = —Q(T m )e\ kT ~) = 

gj 



3h 



7(C 18 0) Q(T ex )e 



8tt 3 Sfi 2 J. (Tex) - MTbg) e h »/ kT e* - 1 rj u 



fl) 



where: 7(C 18 0) is the integrated intensity of the line; Nj 
is the column density of the upper level; gj, Ej and S are 
statistical weight (=2 J + 1), energy of the upper level and 
line strength, respectively; Q(T ex ) is the partition function 
at the temperature T ex ; v the line rest frequency; J^(T ex ) 
and Jv(Tbg) the equivalent Ray leigh- Jeans temperatures at 
frequency v computed for the excitation and background 
temperature (Tbg 00 2.7 K), respectively; \i the molecule's 
dipole moment (0.1102 Debye for C 18 0); r\ v the beam filling 
factor. The total integrated intensity 7(C 18 0) was derived 
from the integral over all the channels with emission instead 
of the area given by the Gaussian fits to take into account 
also non- Gaussian features. 

As for the derivation of the N2H + column density, we 
assumed a unity filling factor for all sources, so that N c i8 Q is 
a beam-averaged value. As excitation temperature, we have 
used the kinetic temperatures listed in Col. 3 of Table 



(i.e. we are assuming LTE conditions for C 18 (3-2)). The 
resulting column densities are listed in Table [9] To check if 
the assumption of optically thin lines is reasonable, we have 
estimated the line optical depths from the solution of the 
line radiative transfer equation: 



-In 1 



Ju(T ex ) — Jv(Tbg) 

and found values smaller than ~ 0.6, consistent with our 
assumption of optically thin lines. 

The CO depletion factor, /d, is defined as the ratio 
between the 'expected' abundance of CO relative to H2 



/d 



and the 'observed' value: 



A ni8n 



(2) 



X° 18Q is the ratio between the observed C 18 column den- 
sity and the observed H2 column density. For this latter, 
we have used the value derived from the 1.2 mm continuum 
peak flux, iVp (H2) (Col. 7 of Table [8]), which is averaged 
over a beam comparable to that of the C 18 observations 
(24" against 19"). 

To compute X^ 18Q , we have taken into account the 
variation of atomic carbon and oxygen abundances with 
distance to the Galactic Cente r follo wing the same proce- 
du re ad opted in Fontani et al. (2006; see also Miettinen et 
al. l201lh . Assuming the canonical abundance of ~ 8.5 x 10 -5 
for the abundance of the main CO isotopologu e in th e neigh- 
bourhood of the s olar system (Frerking et al. Il982k see also 
Langer et al. 1989 and Pineda et al. 2008), we have computed 
the expected CO abundance at the Galactocentric distance 
(Dqc) of each source according to the relationship: 

X^i6 = 8.5 x 10" 5 exp (1.105 - 0.13 £> G c(kpc)) , (3) 

which has been derived according to the abundance gradi- 
ents in the Galactic Disk for 12 C/H and 16 0/H listed in 
Table 1 of Wilson & Matteucci (|l992h , and assuming that 
the Sun has a distance of 8.5 kpc from t he Galactic Centre. 
Then, following Wilson & Rood ([19941 ), we have assumed 
that the Oxygen isotope ratio 16 0/ 18 depends on Dqc ac- 
cording to the relationship 16 0/ 18 0=58.8x D G c (kpc) +37.1, 
which gives: 



xS 18 



C^O 



(58.8L>gc + 37.1) 



(4) 



The results are listed in Table [9] Cols. 2, 3 and 4 
list the C 18 (3-2) integrated line intensity (7(C 18 0)), the 
C 18 column density (N(C 18 0)) and the observed C 18 
abundance (X° 18o ); Cols. 5 and 6 give the expected C 18 
abundance (X^ 18Q ) calculated at the Galactocentric dis- 
tance of the source (Table [1]) and the CO depletion factor 
(/d = Xci8 /X° 18Q ), respectively. 

The mean fr> is ~ 32 and the median value is 29. These 
values are remarkably higher than those obta ined in other 
IRD Cs (P illai et al. l2007l . Miettinen et al. l201lL Hernandez 
et al. l201lh . and are among the highest obtained both in low- 
mass starless cores (Crapsi et al. I2005L Tafalla et al. l2Q06h 
and in massive clumps fr om o bservations with l ow- angular 
resolution (Fontani et al. I2006L Chen et al. l201lh . 

Such a difference could be explained by the fact that 
the transition used in this work to derive fr> has a criti- 
cal density of ~ 5 x 10 4 cm -3 , comparable to that of the 
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Table 8. Parameters derived from the 1.2 mm continuum emission: 1.2 mm continuum flux density (S u ) integrated over 
the 3a rms contour level, clump angular diameter (0 S ), gas mass (M), H2 volume (nn 2 ), column (iV(H2)) and mass surface 
density (S(H2)). The last column lists the mass threshold, Mt h r , fo r a cloud that can form massive stars based on the 
empirical mass-radius relation proposed by Kauffmann & Pillai (2010, see text). 



source 


S v a 

(Jy) 


6 s b 
(") 


M 

(M ) 


7iH 2 

(xlO 4 cm" 3 ) 


iVt(H 2 ) c 
(xlO 22 cm" 2 ) 


N p (R 2 ) d 
(xlO 23 cm" 2 ) 


S(H 2 ) 
(g cm" 2 ) 


M thr 
(M ) 


08477-4359cl 


1.8 


35.6 


86.73 


11.5 


11.0 


1.42 


0.24 


73 


13039-6108c6 


1.04 


40.3 


101.5 


3.90 


5.64 


0.68 


0.12 


127 


13560-6133c2 


0.77 


23.8 


426.8 


6.27 


12.5 


1.09 


0.27 


194 


14183-6050c3 


1.06 


39.6 


207.7 


2.96 


5.95 


0.68 


0.13 


196 


15038-5828cl 


0.49 


24.7 


250.6 


4.62 


8.53 


0.85 


0.19 


175 


15278-5620c2 


1.5 


40.1 


243.7 


3.34 


6.81 


0.82 


0.15 


200 


15470-5419cl 


1.16 


24.2 


310.2 


11.0 


16.4 


1.37 


0.36 


131 


15470-5419c3 


3.06 


54.1 


743.4 


2.37 


7.84 


1.11 


0.17 


382 


15557-5215c2 


2.89 


41.3 


633.4 


3.67 


9.96 


1.55 


0.22 


293 


15557-5215c3 


0.51 


35.8 


194.3 


1.73 


4.07 


0.49 


0.09 


242 


16061-5048cl 


2.1 


28.1 


284.3 


9.54 


14.4 


1.66 


0.31 


134 


16061-5048c4 


1.57 


62.8 


504.2 


1.52 


5.12 


1.22 


0.11 


391 


16093-5128c2 


0.67 


39.6 


481.6 


1.19 


4.29 


0.92 


0.09 


427 


16093-5128c8 


0.86 


40.1 


642 


1.52 


5.57 


0.68 


0.12 


434 


16164-4929c3 


0.35 


36.2 


54.79 


2.28 


3.21 


1.50 


0.07 


122 


16254-4844cl 


1.15 


20.3 


164.2 


17.4 


17.9 


1.43 


0.39 


81 


16435-4515c3 


0.51 


17.7 


147 


30.9 


25.3 


1.20 


0.55 


61 


16482-4443c2 


0.24 


< 24 e 


59.08 


> 3.81 e 


> 4.63 e 


0.66 


> 0.10 e 


< 101 e 


16573-4214c2 


0.93 


7.29 


108.3 


553 


157 


1.89 


3.4 


14 


17040-3959cl 


0.81 


16.6 


3428 


5.80 


23.7 


1.60 


0.52 


505 


17195-3811c2 


2.61 


44.3 


597.9 


5.12 


12.2 


1.77 


0.27 


246 


a from Beltran et al. ([20061): 



b full width half maximum angular diameters deconvolved with the telescope beam, from Beltran et al. (2006); 
c computed from the total clump mass (see text); 
d computed from the 1.2 mm continuum peak flux (see text); 

e unresolved in the SIMBA map. Therefore, tih 2 , iVt(H2) and E(H2) are lower limits, and M t h r is an upper limit. 



central region where the CO freeze-out is expected to be 
important. In fact, the CO freeze-out timescale, which de- 
pends on the H2 volume density, becomes shorter than the 
free-fall timescale at densities larger th an a few 10 4 cm -3 
(see e.g. Bergin & Tafalla 120071 Caselli l201lh . On the con- 
trary, in most of the studies performed so far /d was derived 
from the CO (1-0) or (2-1) lines, which trace lower-density 
gas where CO is significantly non-depleted, resulting in val- 
ues o f ,fp s maller than o r com parable to ^ 10 ( e.g. Crapsi 
et al. 120051 . Tafalla et al. l2006l . Pillai et al. l2007h . However, 
the discussion of this result requires three main comments. 
First, the values of the 'canonical' CO abundance measured 
by other authors in different objects vary by a factor of 2 
(see e.g. Lacy et al. 1994; Alves et al. 1999), but because 
the value assumed by us was d erived from the Taurus star 
forming regions (Frerki ng et al. 1 19821 ) and used in previous 
estimates (Crapsi et al. 120051 : Emprechtinger et al. 1 20091 ), it 
is likely the most appropriate to make comparison with low- 
mass dense cores. Second, the angular resolution of our ob- 
servations allows us to derive only average values of fr> over 
angular regions much larger than the typical fragmentation 
scales (few arcseconds at the distance of our targets) , which 
in reality may have complex structures, so that our estimates 
should be considered as lower limits. Finally, high depletion 
of C 18 could be due to mechanisms other than freeze-out 
in regions where young stellar objects are embedded, and 
can decrease the CO abundance even in warm e nviron ments 
(for a detailed explanation, see Fuente et al. l2012h . This 



comment especially concerns the clumps where the average 
kinetic temperature exceeds the CO sublimation tempera- 
ture of - 20 K. 



5 DISCUSSION 

5.1 Virial state of the clumps 

In order to investigate if our targets are gravitationally 
stable, we have computed the virial masses of the clumps 
assuming virial equilibrium and negligible contributions of 
magnetic field and surface pressure. Assuming also that the 
cores are spherical, one can demonstrate that the gas mass 
is given by: 

M V ir(M ) ~ kAv 2 (km s- 1 )i?(pc) , (5) 

where R is the clump radius, Av 2 the line width and k is a 
multiplicative factor that depends on the gas density distri- 
bution as a function of the dista nce to the clump centre (see 
Eq. (3) in MacLaren et al. 1988). For a homogeneous clump 
(i.e. p = const.), k ~ 210, while for p ~ r -2 , k ~ 126. The 
millimetre continuum maps do not allow to derive the den- 
sity profile of the clumps because of the low angular resolu- 
tion and sensitivity to extended emission. However, previous 
studies suggest that the density structure of massive clumps 
reasonably can be approximated by a constant density in 
the inner regions, and by a steep power-law in the outer lay- 
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Table 9. Parameters used to determine the CO depletion factor: integrated intensity of the C ls O (3-2) 
line (7(C 18 0))), C ls O total column density (7V(C ls O)), observed C ls O abundance (X° 18q ), 'expected' 
C ls O abundance (X^ 18Q ), and CO depletion factor (/d = ^qiso/^c 18 0^' ^ e errors are gi ven between 
parentheses. The uncertainty on /d based on the calibration errors affecting the C ls O and H2 column 
densities is of the order of the 40-50%. 



source 


7(C 18 0) 
(K km s- 1 ) 


N(C ls O) 
(xl0 15 cm- 2 ) 


X° 
c 18 o 

(xlO- 8 ) 


c 18 o 
(xlO- 7 ) 


/d 


08477-4359cl 


6.5(0.1) 


3.03(0.05) 


2.14 


1.45 


7 


13039-6108c6 


1.5(0.2) 


0.65(0.09) 


0.95 


2.08 


22 


13560-6133c2 


1.8(0.1) 


1.08(0.06) 


0.99 


2.68 


27 


14183-6050c3 


0.9(0.1) 


0.41(0.04) 


0.60 


2.54 


42 


15038-5828cl 


0.89(0.09) 


0.52(0.05) 


0.61 


3.30 


54 


15278-5620c2 


2.5(0.1) 


1.10(0.04) 


1.34 


2.93 


22 


15470-5419cl 


2.3(0.1) 


1.36(0.06) 


0.99 


3.46 


35 


15470-5419c3 


2.6(0.1) 


1.06(0.04) 


0.96 


3.46 


36 


15557-5215c2 


4.2(0.1) 


1.90(0.04) 


1.20 


3.79 


32 


15557-5215c3 


1.6(0.1) 


0.76(0.05) 


1.57 


3.79 


24 


16061-5048cl 


8.8(0.2) 


4.7(0.1) 


2.83 


3.39 


12 


16061-5048c4 


3.0(0.1) 


1.23(0.04) 


1.01 


3.39 


34 


16093-5128c2-a 


3.55(0.08) 


1.59(0.04) 


1.72 


5.09 


29 


16093-5128c8 


1.00(0.06) 


0.44(0.03) 


0.65 


5.09 


78 


16164-4929c3 


3.01(0.07) 


1.40(0.03) 


0.93 


2.80 


30 


16254-4844cl 


2.82(0.09) 


1.93(0.06) 


1.35 


3.39 


25 


16435-4515c3-a 


1.73(0.09) 


0.68(0.07) 


0.57 


4.13 


73 


16482-4443c2 


4.7(0.1) 


2.96(0.06) 


4.45 


3.87 


9 


16573-4214c2 


7.8(0.1) 


2.2(0.3) 


1.18 


2.98 


25 


17040-3959cl 


6.32(0.09) 


5.32(0.08) 


3.32 


1.59 


5 


17195-3811c2 


3.07(0.08) 


1.33(0.03) 


0.75 


4.04 


54 



In Fig. [3] we compare the mass derived from the dust 
continuum emission to Mvir obtained assuming p ~ const. 
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Figure 3. Virial masses against gas masses computed from the 
1.2 mm dust continuum emission. Mvir is calculated assuming a den- 
sity distribution of the type p = const, (left panel) and p ~ r~ 2 (right 
panel). In both panels, the line indicates Mvir= M. 



(left panel) and p 



(right panel). We note that Mvir 



is on average larger than M for the case p ~ const., while it 
is on average smaller for the other case. Because the overall 
density distribution of the clumps is likely in between these 
two extreme cases, we suggest that, on average, the clumps 
are close to virial equilibium. Our findings are in good agree- 
ment with those derived by Lopez-Sepulcre et al. (|2010h 
in a sample of high-mass clumps supposed to span a wide 
range of evolutionary stages. The gas masses of the high- 
mass clumps studied by Lopez-Sepulcre et al. (|2010|) were 
derived from millimetre continuum measurements, and the 
dust temperatures assumed are in agreement with ours (see 
their Table 1). Moreover, they estimated the virial masses 
from C 18 (2-1) line widths and assuming p ~ const. As 
shown in Fig. [4l we do not highlight systematic differences 
between the two sub-samples, and between infrared-dark 
and infrared-bright sources either. On average, the sources 
of the Lopez-Sepulcre et al.'s sample appear to have slightly 
larger virial masses, so that potentially they might be less 
gravitationally bound than the clumps studied in this work. 
However, we stress that the parameters from which the mass 
estimates are derived are affected by large errors (especially 
the dust mass opacity, the distance and the gas-to-dust ra- 
tio) and difficult to quantify. For example, the dust mass 
opacity coefficient ca n intr oduce a factor 2 in the uncer- 
tainty (Beuther et al. 2002). Hence our conclusions on the 
virial stability of the targets must be taken with caution. 



ers (see e.g. Beuther et al. 12001 Fontani et al. l200l Hill et 
al. l2005h . 
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p = const 



® 3 



a □ a, 



o 



tend to have lower mass (average M of ^ 260 M© from 
the 1.2 mm continuum), smaller source- averaged H2 column 
density (average Nt (H.2) of 8.9 x 10 22 cm -2 ), and smaller 
H2 surface density (average £(£[2) of 0.19 gr cm -2 ) than 
the infrared-bright clumps (average M of ^ 370 M©, av- 
erage iVt(H 2 ) of 1.3 x 10 23 cm" 2 , average E(H 2 ) of 0.27 g 
cm" 2 ). 

The fact that on average the group of the 24 /im-dark 
sources has lower Tk and higher CO depletion factor than the 
24 /im-bright objects suggests that it likely contains objects 
less evolved, and this is in accordance with other compara- 
tive studies of massive clumps with or with out indications 
of embedded star formation (Hill et al. 120051 ). However, the 
dark sources also have lower M and H2 column and surface 
density, and this could indicate that the 24 pm dark clumps 
may be destined to form less massive stars/clusters, thus 
explaining the non-detection at mid-infrared wavelengths, 
although material could still be accreting. 



12 3 4 

Log 10 [M (M )] 

Figure 4. Virial masses against gas masses computed from the dust 
continuum emission for the clumps studied in this work (tria ngles) and 
the massive clumps observed by Lopez-Sepulcre et al. (|201Cl squares). 
Filled and open triangles represent sources detected and undetected 
at 24 /im, respectively, while the filled and open squares represent the 
infrare d-brig ht and infrared-dark sources selected by Lopez-Sepulcre 
et al. |2010), respectively. The cross indicates the clump with no 24 
/im image available, 13039-6108c6. Mvir is calculated assuming ho- 
mogeneous clumps. The line indicates Mvir= M. 

5.2 24 pm dark versus 24 /im bright clumps 

The presence of mid-infrared emission in molecular clumps 
is very often the sign of embedded star formation activity. 
Therefore, the IRDCs detected at 24 /im could harbour ob- 
jects more evolved than those undetected. If this is the case, 
the two groups should have physical properties indicative 
of a different evolutionary stage. In Fig. [5] we show his- 
tograms which compare some physical and chemical prop- 
erties of the IRDCs detected and undetected at 24 /im: line 
widths, kinetic temperature, column (N t (H2)), volume and 
surface density of H2, gas masses (both from continuum and 
the virial theorem), CO depletion factor and clump diam- 
eter. For Mvir, we have considered the values calculated 
assuming p ~ const, bearing in mind that those obtained 
in the case p ~ r~ 2 are just systematically lower by a fac- 
tor ~ 1.7. For the unresolved source 16482-4443c2, we have 
not included the n(H2) and S(H2) lower limits, and replaced 
Nt (H2) with Np(H2). Although the statistics is poor because 
the two sub-samples contain eight 24 /im bright and twelve 
24 /im dark objects, the comparative histograms are useful 
to check if one can notice clear systematic differences. An 
inspection of the comparative histograms indicates that the 
24 /im dark clumps have lower Tk and higher /d- In fact, 
the average Tk and /d of the infrared-bright group turns 
out to be ~ 20 K and 28, respectively, while that of the 
infrared-dark group is ~ 16 K and 35, respectively. Also, 
if we exclude from the statistical analysis the clear outliers 
17040-3959cl and 16573-4214c2, the 24 pm dark clumps 



5.3 Relation between Tk and line widths 

Two indicators of active star formation are the line width 
and the gas temperature, which are both expected to be- 
come higher with increasing star formation activity. In fact, 
the line widths of dense gas tracers are found to be higher 
in more evolved star-forming clumps than in quies cent re- 
gions of infrared-dark clouds (e.g. Hill et al. 12010 . Ragan 
et al. 120121 ). while the cores associated with protostars are 
on average warmer than the starless cores, both in low- and 
high-mass star for ming regions (e.g. F oster e t al. [2009, Em- 
prechtinger et al. 2009|, Ragan et al. 1201 lh . Fig. [6] shows 
that there is a slight correlation between Tk and the C 18 
line width. If we exclude from the statistical analysis 15038- 
5828cl, for which Ay is much higher than that of the other 
targets and has a spectrum with poor S/N (see Fig. IB2|) . a 
linear least square fit to the data yields a slope of ~ 0.09. 
Statistical correlation between Tk and line width can be in- 
vestigated also through non-parametric statistical methods, 
like the Kendall's r correlation coefficient. This measures 
the rank correlation, namely how two quantities are ordered 
similarly when ranked by each of therrjf|. r can range be- 
tween 1 (perfect agreement between the two rankings) and 
— 1 (one ranking is the reverse of the other). We find r ~ 0.23 
between Tk and line width. These results suggest a faint cor- 
relation between the two parameters, which indicates that 
the warmer clumps tend to also have larger line widths, i.e. 
tend to be more turbulent. However, the large dispersion 
cannot allow us to draw any firm conclusion. 



5.4 Relation between CO depletion and other 
physical parameters 

One of the main results of this work is the high CO deple- 
tion factor measured in our targets, with values higher than 
those foun d in comparable high-mass clumps (e.g. Mietti- 
nen et al. l201lL Chen et al. l201lh , a nd ev en in low-mass 
pre-stellar cores (see e.g. Crapsi et al. l2005h . The fact that 
we have derived /d from the C 18 (3-2) transition, which 
has a critical density of 5 x 10 4 cm -3 , certainly makes the 



for details see, e.g., http:/ /www. statsoft.com/textbook/nonparametric-statistics/ ; 
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Figure 5. Histograms comparing some of the physical and chemical properties of the clumps detected (solid line) and 
undetected (dashed line) in the Spitzer-MIPS 24 /im image (see Figs. IAl|) . In the panels showing iV(H2), n(H2) and £(H2) 
we have not included the lower limits calculated for 16482-4443c2. 



difference with respect to previous works in which /d was 
calculated from the lower excitation transitions tracing ma- 
terial where freeze-out of CO is expected to be less impor- 
tant. In fact, freeze-out of CO and other neutrals onto dust 
grains is favoured in gas characterised by low temperatures 
(T ^ 20 K) and high densities (n ^ 10 5 cm -3 ), in which 
the freeze-out timescale is generally shorter than the free- 
fall timescale (e.g. Bergin & Tafalla l2007h . In fact, high CO 
depletion factors (of the order of 10 or more) were measured 
towards the dense and cold nucle i of low-mass pre-stellar 
cores (see e.g. Crapsi et al. 2005). After the formation of 
the protostellar object (s) at the core centre, evaporation of 
C O start s and the CO depletion factor drops (Caselli et 
al. [2002b). This theoretical prediction has been partly con- 
firmed by observations of both low- and high-mas s den se 
cores (e.g. Emprechtinger et al. 120091 Fontani et al. 1 20061 ). 



In Fig. [3 we show fr> against Tk (in the left panel) and 
n(H2) (in the right panel). At a first glance, the two plots do 
not reveal clear (anti-) correlations between the parameters. 
As made in Sect. 15.31 we have performed a closer inspection 
of the data using the non-parametric ranking statistical test 
Kendall's r (see Sect. 15. 3|) . If we consider all points, we find 
a faint anti-correlation between /d and Tk (Kendall's r — 
-0.34) and also between /d and n(H2) (Kendall's r — - 
0.12). The anti-correlation between Tk and /d is consistent 
with the increase of CO depletion with decreasing temper- 
ature, as already found in both low- and hig h-mas s dense 
cores (Emprechtinger et al.|2009|, Fontani et al. |2006T L while 
that between fr> and n(H2) is the opposite of what expected 
from chemical models. Even if we exclude from the statisti- 
cal analysis 16573-4214c2, which has n(H2) markedly much 
higher than that of the other members of the sample, the 
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Figure 6. Gas kinetic temperature, Tk, derived from ammonia, 
against the C ls O line widths, AV. The straight line is a least square fit 
to the data excluding 15038-5828cl, the clump with AV > 5 km s _1 , 
i.e. much higher than any other clump and based on a very noisy spec- 
trum (see Fig, IB2|) . The slope of the linear fit is ~ 0.09. Filled and open 
symbols represent clumps detected and undetected at 24 ^m, respec- 
tively. The cross indicates the clump with no 24 fim image available, 
13039-6108c6. Typical errorbars are depicted in the top-right corner. 



correlation between fr> and n(H2) is not observed (Kendall's 
r = -0.11). We point out, however, that all the parame- 
ters obtained are average values over angular regions much 
larger than the size expected for the embedded condensa- 
tions, while the correlations are predicted for single cores. 
Moreover, clumps embedded in different clouds may be af- 
fected by different environmental conditions, so increasing 
the dispersion. In fact, Emprechtinger et al. (2009) found 
that dense cores in Perseus showed the best trends between 
the physical parameters, whereas when including cores from 
other star-forming regions the dispersion was significantly 
larger. 

Finally, Fig. [8] suggests that /d is not clearly (anti- 
) correlated to either the clump mass M and and the ra- 
tio Mvir/M which should show the tendency of a clump 
to be stable against gravitational collapse. When excluding 
the clear outlier 17040-3959cl (whose mass is statistically 
much larger than the rest of the sample) , a slight correlation 
is found between /d and M (Kendall's r = 0.4). However, 
again the large dispersion of the data and the big errors on 
both parameters do not allow us to claim firm conclusions. 



6 CONCLUSIONS 

We have undertaken a molecular line study, through obser- 
vations of C 18 (3-2), NH 3 (1,1) and (2,2), N 2 H+ (3-2) 
and N 2 D + (4-3), of 21 IRDCs in the southern hemisphere 
with the aim of characterising the initial conditions of the 
high-mass star and cluster formation process. The sample 
targeted was selected from the high-mass millimetre clumps 
not detected in the MSX images identified by Beltran et 
al. (2006), and includes sources with and without Spitzer 
24 fim emission. The NH3 inversion transitions have been 
observed with the ATCA. The rotational transitions of the 
other molecular species have been observed with the APEX 
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Figure 8. CO depletion factor (/d) against gas mass 
derived from the 1.2 mm continuum (top panel) and 
Mvir/M mass ratio (bottom panel). The symbols 
have the same meaning as in Fig. [6] 

Telescope. We have detected C 18 and ammonia emission in 
all clumps, and N 2 H + emission in all the 12 sources observed 
in this line. Only one source has been marginally detected in 
N 2 D + (4-3), which is, to our knowledge, the first detection 
of this line in an IRDC. The clumps have a median mass of 
~ 244 Mq, appear to be gravitationally bound and possess 
mass, H2 column and surface densities consistent with being 
potentially the birthplace of high-mass stars. The most strik- 
ing result of the work is the high average CO depletion factor 
(derived from the expected C 18 0-to-H2 column density ra- 
tio compared to the observed value), which is in between 5 
and 78, with a mean value of 32 and a median of 29. These 
values, derived from the C 18 (3-2) transition (which traces 
gas denser than the lower- excitation lines commonly used in 
previous works), are larger than the typical CO depletion 
factors measured towards other IRDCs, and are comparable 
to or larger than the values derived in the low-mass pre- 
stellar cores closest to the onset of gravitational collapse. 
Also, a faint ant i- correlation is found between /d and the 
gas kinetic temperature. These results suggest that the ear- 
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Figure 7. CO depletion factor (/d) against gas kinetic temperature (left panel) and H2 volume density (right panel). The 
symbols have the same meaning as in Fig. [6] 



liest phases of the high-mass star and stellar cluster forma- 
tion process are characterised by CO depletions larger than 
in low-mass pre-stellar cores. The clumps have an average 
temperature around 17 K. We have found marginal statisti- 
cal differences between the clumps detected and undetected 
in the 24 /im Spitzer images, with the Spitzer-dark clumps 
being on average colder, less massive and having lower H2 
column and surface densities, but higher CO depletion fac- 
tors. This indicates that the Spitzer-dark clumps are either 
less evolved or destined to form stars and stellar clusters less 
massive than the Spitzer-bright ones. 
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APPENDIX A: SPITZER AND SIMBA IMAGES 
APPENDIX B: SPECTRA 

This paper has been typeset from a TgX/ LT^X file prepared 
by the author. 
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Figure Al. Superimposition of the Spitzer-MIPS 24 /im images on the SIMBA 1.2 mm continuum maps (blue contours: 
first contour and step correspond to the 3 a rms level) towards the IRDCs studied in this work and observed with Spitzer- 
MIPS (for clump 13039-6108c6, the Spitzer-MIPS image is not available). In each frame, the red circle indicates the targeted 
millimetre clump. 
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Figure Bl. Spectra of the twelve sources observed and detected in both N2H+ (3-2) and C ls O (3-2). In each frame, we 
show a range of ±20 km s _1 around the systemic velocity adopted to centre the spectra (given in Table [T} and derived from 
CS observations. For two objects, 16093-5128c2 and 16435-4515c3, the spectrum contains velocity components separated 
by several km s _1 , so that we show a broader velocity range. 
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Figure B2. Spectra of the nine sources observed and detected in C ls O (3-2) only. As in Fig. IB1I the velocity interval on 
the x-axis range from Vlsr ~~ 20 km s _1 to Vlsr + 20 km s _1 . 



